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ABSTRACT 

Context. The study of detached eclipsing binaries which are members of stellar clusters is a powerful way of determining the properties of the 
cluster and of constraining the physical ingredients of theoretical stellar evolutionary models. 

Aims. DWCarinae is a close but detached early B-type eclipsing binary in the young open cluster Collinder228. We have measured accurate 
physical properties of the components of DW Car (masses and radii to 1%, effective temperatures to 0.02 dex) and used these to derive the age, 
metallicity and distance of Collinder228. 

Methods. The rotational velocities of both components of DW Car are high, so we have investigated the performance of double-Gaussian 
fitting, one- and two-dimensional cross-correlation and spectral disentangling for deriving spectroscopic radial velocites in the presence of 
strong line blending. Gaussian and cross-correlation analyses require substantial corrections for the effects of line blending, which are only 
partially successful for cross-correlation. Spectral disentangling is to be preferred because it does not assume anything about the shapes of 
spectral lines, and is not significantly affected by blending. However, it suffers from a proliferation of local minima in the least-squares fit. We 
show that the most reliable radial velocities are obtained using spectral disentangling constrained by the results of Gaussian fitting. Complete 
Stromgren uvby light curves have been obtained and accurate radii have been measured from them by modelling the light curves using the 
Wilson-Devinney program. This procedure also suffers from the presence of many local minima in parameter space, so we have constrained the 
solution using an accurate spectroscopic light ratio. The effective temperatures and reddening of the system have been found from Stromgren 
photometric calibrations. 

Results. The mass and radius of DWCarA are Ma = 11.34 ± 0.12Me and Rx = 4.558 ± 0.045Ro. The values for DWCarB are 
Mb = 10.63 ± 0.14 Mo andi?B = 4.297 + 0.055 Rq. Stromgren photometric calibrations give effective temperatures of T^s/^ = 27 900 + 1000 K 
and Teffg = 26 500 ± 1000 K, and a reddening of Et-y = 0. 18 + 0.02, where the quoted uncertainties include a contribution from the intrinsic 
uncertainty of the calibrations. The membership of DWCar in Cr228 allows us to measure the distance, age and chemical composition of the 
cluster. We have used empirical bolometric corrections to calculate a distance modulus of 12.24 ± 0.12 mag for DWCar, which is in agreement 
with, and more accurate than, literature values. A comparison between the properties of DWCar and the predictions of recent theoretical 
evolutionary models is undertaken in the mass-radius and mass-Feff diagrams. The model predictions match the measured properties of 
DWCar for an age of about 6Myr and a fractional metal abundance of Z a: 0.01. 
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- stars: emission-line. Be - stars: distances - stars: individual: DW Carinae - open clusters and associations: individual: Collinder 228 



1. Introduction 

The study of detached eclipsing binary star systems (dEBs) 
is one of the few ways in which the absolute properties of 
stars can be measured directly and accurately (Andersen 1991> . 
These objects are therefore of fundamental importance to the 

Send offprint requests to: J. Southworth 

* Based on observations carried out with the Stromgren Automatic 
Telescope (SAT) and the FEROS spectrograph at ESO, La Silla, Chile 
(62.H-0319, 62.L-0284, 63.H-0080, 64.L-0031) 



understanding of the characteristics of single stars. The masses 
and radii of the component stars in dEBs can be determined to 
accuracies of better than 1 % from the analysis of light curves 
and double-lined radial velocity curves (e.g., Southworth et al. 
I2005bl Torres et al. 1997 1. The effective temperatures of the 
two stars can be found by spectral analysis or photometric cal- 
ibrations, allowing the luminosities and distance of the stars 
to be calculated directly (e.g., Southworth, Maxted & Smalley 
12005a. Clausen & Gimenez .l991< . 
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One of the most important uses of the physical properties 
of dEBs is as a check of the predictions of theoretical mod- 
els of stellar evolution. The two components of a dEB have 
the same age and chemical composition, as they were born 
together, but in general different masses, radii and luminosi- 
ties. Theoretical models must therefore be able to match their 
accurately-known properties for one age and chemical com- 
position. A good match is often achieved (e.g., Southworth, 
Maxted & Smalley"2004b'i for two reasons. Firstly, the models 
are in general fairly reliable, and the effects of improved in- 
put physics lead to only minor adjustments in the evolutionary 
phases for which we are able to study dEBs. Secondly, there 
are several quantities which are not in general independently 
known for dEBs, including age, metal abundance and helium 
abundance, and these can be freely adjusted to help the models 
match the observed properties of a dEB. This last point means 
that it is extremely difficult to investigate the success of the 
treatment in theoretical models of several physical phenomena, 
for example convective core overshooting, mixing length, mass 
loss and opacity (see Cassisi 120041 120051 . 

Extra constraints on the properties of a dEB can be 
obtained if it is a member of a co-evolutionary group of 
stars (Southworth, Maxted & Smalley 2004a, Hebb, Wyse & 
Gilmore r2004> . If the age and chemical composition of the stars 
is known independently then these quantities can be used to 
constrain the models, making it more challenging for them to 
match observations and so allowing the effects of more subtle 
physical phenomena to be investigated (e.g., Torres & Ribas 
I2002> . If the age and composition of the stellar group is not 
known then these may be found very precisely by comparing 
model predictions to the properties of the dEB (Southworth et 
al.|2004a 2004b Thompson et al. 2001 1. In addition, dEBs are 
excellent distance indicators (e.g., Guinan et al. ll998l Clausen 
12004 1, allowing accurate and precise distances to be found 
to the parent stellar systems in several empirical and semi- 
empirical ways (Southworth et al. l2005al Guinan et al. ll998> . 

1.1. DW Carinae 

Here we present a study of the young, early B-type system 
DWCar, which is a member of the Southern open cluster 
Collinder228 (Ferrer et al. 1985). We have obtained complete 
and accurate Stromgren uvby light curves, which are presented 
in Clausen et al. (;2006 hereafter Paper I) along with a re- 
vised orbital ephemeris. We have also obtained extensive spec- 
troscopy using the FEROS echelle spectrograph and the 1.5 m 
telescope at ESO La Silla. 

DW Car (TableQ is composed of two similar B 1 V stars in 
a 1.3 day orbit. The stars are well separated, despite the short- 
ness of the orbital period, as they are essentially unevolved. 
Such a system is particularly useful because there are very few 
well-studied components of dEBs with masses of IOMq or 
greater which are close to the zero-age main sequence (ZAMS) 
(Andersen 1 199 11 Gies|2003). High-mass unevolved dEBs can 
be particularly useful for investigating the effects of the opaci- 
ties used in theoretical models (Torres et al. .l997j . 



Table 1. Identifications and astrophysical data for DWCar. To 
is a reference time of mid-eclipse of the primary star. 
References: (1) Cannon n^25^; (2) H0g et al. (^2000^; (3) 
Levato & Malaroda ( .1981 K (4) Paper I; quantities in parenthe- 
ses are uncertainties in the final digit of the quanities. 





DWCar 


References 


Henry Draper Catalogue 


HDE 305543 


1 


Tycho number 


TYC 8957-1314-1 


2 


C2000 


10 43 10.1 


2 


(52000 


-60 02 12 


2 


Spectral type 


BlV-l-Bl V 


3 


V 


9.675 ± 0.005 


4 


b - y 


0.067 ± 0.005 


4 


mi 


0.046 ± 0.008 


4 




-0.017 ± 0.008 


4 


P 


2.528 ± 0.007 


4 


To (HJD) 


2446 828.6692(4) 


4 


Orbital period (days) 


1.32774934(44) 


4 



The eclipsing nature of DWCar was discovered by 
Hertzsprung (1924i, who measured the orbital period to be 
0.66382 d from a photographic light curve containing 251 ob- 
servations, van den Hoven van Genderen ('1939^ refined the 
orbital period using a 914-point photographic light curve, but 
stated that the period should be doubled even though there was 
no clear evidence that adjacent minima had different depths. 
Gaposchkin ( 1952 1953 1 published a photographic light curve 
containing over one thousand points and found an orbital pe- 
riod of 1.3277504d. His estimated masses and radii are quite 
close to the values we measure below. 

A spectroscopic orbit for DW Car was published by Ferrer 
et al. (1985 1 based on 67 photographic observations. They 
found no evidence for orbital eccentricity, and a systemic ve- 
locity which is consistent with membership of Cr228. A spec- 
tral classification ofBlV + BIV was provided by Levato & 
Malaroda (1981 1. 

We shall refer to the primary and secondary components 
of DW Car as star A and star B, respectively, where star A is 
eclipsed by star B at phase 0.0. Star A has a greater surface 
brightness and mass than star B. 

1.2. Collinder228 

The open cluster Cr 228 is a young and sparse cluster located 
in the Carina spiral arm feature. Its immediate surroundings are 
quite nebulous, making photometric studies very difficult. 

Feinstein et al. ( 1976 1 found a reddening of Eb-v - 0.33 + 
0.08 mag, an age of less than 5 Myr and a distance modulus of 
Vo - My - 12.0 + 0.2 mag on the basis of a photoelectric UBV 
study. 

The, Bakker & Antalova |l980l obtained photoelectric ob- 
servations in the Walraven system and found a distance modu- 
lus of 12.0 mag assuming a normal reddening law (or 1 1 .7 mag 
assuming Ry - - 4-0). 

Turner & Moffat ( 1980 1 obtained photoelectric UBV obser- 
vations and showed that the reddening in the region was normal 
and that the open clusters Trumpler 14, Trumpler 15 and Cr228 
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are at a common distance modulus of 12.15 + 0.14mag. They 
found that the age of Cr 228 must be at least 1 Myr from the 
presence of main sequence stars at {B - V)o = -0.1. 

Tapia et al. (1988) obtained JHK photometry of the re- 
gion, and found the reddening characteristics of Cr228, and 
other nearby clusters, to be anomalous. These authors found 
that the absolute visual extinction. Ay, is quite different when 
calculated from the colour excess Eg-v rather than from colour 
excesses for longer-wavelength passbands. Under the assump- 
tion that Eb-v is anomalous in this region, they found Ay = 
1.94 + 0.16 mag from Ey-K, giving a distance modulus of 
Vo-Mv = 11.6 + 0.4. 

Carraro & Patat J2001t obtained UBVRI photometry of 
Cr228 and several nearby clusters. They found an extinction 
of Eb-v - 0.30 + 0.05 mag, in agreement with Feinstein et al. 
J197 6l but disagreeing with Tapia et al. ( ,1988j . Their distance 
modulus, 11.4 + 0.2 mag, is smaller than that derived in most 
previous studies. Carraro & Patat suggested that about 30% of 
the stars in Cr 228 are binary, in good agreement with the spec- 
troscopic study of Levato et al. ( 1990). 

Massey, DeGioia-Eastwood & Waterhouse ( 2001) obtained 
new spectral types for 16 stars, and collected literature val- 
ues for another 40 stars, in the region of Cr228. They found 
a median Eb-v of 0.37 and a distance modulus of 12.5 mag. 
They obtained the ages of the higher-mass unevolved stars from 
interpolation in effective temperature and absolute bolometric 
magnitude using the theoretical evolutionary tracks of Schaller 
et al. ( 1992). The age and standard deviation is 2.2 + 1.0 Myr, 
but these authors note that the appearance of the HR diagram 
suggests a spread in ages in Cr228. 

It is clear from the above discussion that the physical prop- 
erties of the open cluster Cr 228 are rather uncertain, because 
the cluster is quite sparse and located in a nebulous area. 
The brightest member of this cluster, QZ Car (HD 9326), is a 
quadruple system which contains two spectroscopic binaries, 
one of which shows eclipses of a semi-detached nature (Mayer 
et al.EOOD- 



Table 2. Logbook of the spectroscopic observations of 
DWCar HJD refers to the exposure midpoint, ID to the 
identification number, SD to if it was included in the spec- 
tral disentangling analysis and fexp to the exposure time. The 
column headed 'Obs' indicates the observer (H indicates 
Heidelberg/Copenhagen guaranteed time observations). The 
signal to noise ratio, S/N, of single exposures was measured 
at approximately 4000 and 5000 A. 



HJD minus 
2400000 


ID 


SD 


Obs 


Orbital 

phase 


^exp 

(s) 


S/N 


51144.83060 


si 


y 


H 


0.73512 


600 


35-60 


51145.84400 


s2 




H 


0.49836 


600 


30-55 


51146.86755 


s3 


y 


H 


0.26925 


600 


40-75 


51147.86649 


s4 




H 


0.02161 


600 


35-60 


51148.80902 


s5 


y 


H 


0.73148 


600 


40-85 


51149.87170 


s6 




H 


0.53184 


420 


25^0 


51150.83746 


s7 


y 


H 


0.25920 


600 


40-60 


51152.73657 


s8 




H 


0.68953 


600 


25-60 


51172.69544 


s9 


y 


H 


0.72164 


900 


50-80 


51174.78088 


slO 


y 


H 


0.29229 


900 


75-110 


51176.80734 


sU 


y 


H 


0.81853 


600 


65-95 


51178.84953 


sl2 


y 


H 


0.35661 


600 


70-90 


51179.87450 


sl3 


y 


H 


0.12857 


600 


65-100 


51180.81349 


sl4 


y 


H 


0.83578 


600 


60-80 


51181.87490 


sl5 


y 


H 


0.63519 


600 


65-90 


51182.86220 


sl6 


y 


H 


0.37877 


600 


60-100 


51184.86864 


sl7 




H 


0.88993 


600 


55-75 


51187.87364 


sl8 




H 


0.15316 


600 


35-65 


51189.86682 


sl9 


y 


H 


0.65433 


600 


50-75 


51190.84867 


s20 




H 


0.39382 


600 


50-70 


51191.86279 


s21 


y 


H 


0.15760 


600 


50-75 


51192.80846 


s22 


y 


H 


0.86984 


600 


65-85 


51193.87953 


s23 




H 


0.67652 


720 


65-95 


51197.85496 


s24 


y 


H 


0.67063 


600 


40-65 


51199.87004 


s25 


y 


H 


0.18830 


900 


60-85 


51209.72395 


s26 




JVC 


0.60981 


1800 


110-140 


51382.51044 


s27 


y 


H 


0.74468 


1800 


50-125 


51386.50761 


s28 




H 


0.75517 


1800 


70-135 


51563.86655 


s29 


y 


JVC 


0.33380 


2400 


105-150 



2. Observations and data reduction 

Light curves of DWCar in the Stromgren uvby photomet- 
ric system were obtained from the Stromgren Automated 
Telescope at ESO La Silla. These are presented in Paper I 
and contain 518 observations in each passband. Paper I also 
presents an updated orbital ephemeris and standard Stromgren 
uvby and Crawford Hj3 indices for DW Car. 

2.1. Spectroscopic data 

29 high-resolution spectra were obtained using the FEROS 
fibre-fed echelle spectrograph on the ESO 1.52 m tele- 
scope at La Silla (Kaufer et al. [T9991 ^000) between 1998 
November and 2000 January. This spectrograph is located in 
a temperature-controlled room and covers the spectral region 
from the Balmer jump to 8700 A without interruption and 
at a constant velocity resolution of 2.7kms"'px"' {A/AA = 
48000). An observing log is given in Table|3 



A modified version of the MIDAS FEROS package, pre- 
pared by H. Hensberge, was used for the basic reduction'. This 
package pays careful attention to background removal, defini- 
tion and extraction of the individual orders and wavelength cal- 
ibration, and provides a significant improvement on the stan- 
dard quick-look reduction pipeline. The observations were re- 
duced night by night using calibration exposures (thorium- 
argon and flat field) obtained during the preceding afternoon 
and following morning for each night. Standard (rather than 
optimal) extraction was applied, and no order merging was 
attempted. A standard error of the wavelength calibration of 
about 0.002-0.003 A was typically obtained. 



' see these websites for further details: 
http : //www . Is . eso . org/lasilla/sciops/2p2/E2p211/FEROS/DRS 
http ; //www . Is . eso . org/lasilla/sciops/2p2/Elp5M/FEROS/ Reports| 
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6550 6560 

Wavelength (A) 

Fig. 1. Continuum-normalised Ha profiles of DW Car at phases 
0.29, 0.50 and 0.72 (spectra s 10, s2 and s9). Spectra s2 and 
s 9 have been shifted by +1 and +2, respectively, for clarity. 



3. Radial velocity analysis 

The FEROS spectra of DW Car cover essentially the entire op- 
tical wavelength range, but are disappointingly short of stellar 
features. The only absorption lines which are easy to discern 
are the hydrogen Balmer lines and several He I lines. The high 
rotational velocity of both stars, approximately 170 km s"', 
twinned with the fact that the spectrum of each star is diluted 
by the flux of the companion, means that only the strongest 
spectral lines are visible. 

The Ha line shows significant emission with an absorption 
superimposed on the line core. The emission and emission ab- 
sorption do not follow the orbital motion of the stars but remain 
at the systemic velocity of the system. Slight emission is also 
visible in the centre of H/3 for some of the spectra, depending 
on signal to noise ratio. Fig.^shows Ha spectra at phases 0.27, 
0.50 and 0.72. Whilst the emission appears to be twice as strong 
at phase 0.50, this is simply an artefact of the continuum nor- 
malisation as the system is about half as bright at eclipse mid- 
point compared to during quadrature phases. Further investiga- 
tion of the emission line strengths will require flux-calibrated 
spectra (D. Lennon, private communication). 

Narrow regions around the five He I absorption lines at 
4026, 4388, 4471, 4922 and 6678 A were chosen for our ra- 
dial velocity analysis. The remaining visible He I lines are all 



4920 4930 4940 

Wavelength (A) 



s 27 



6660 6670 6680 6690 6700 

Wavelength (A) 

Fig. 2. Continuum-normalised example spectra in the regions 
of the five He I lines used in the radial velocity analysis. 
Spectrum s 29 (phase 0.33) is shown in the top four panels and 
spectrum s 27 (phase 0.74) in the last panel. Shallow lines can 
also be seen at 4009 A (He I) and 4481 (Mgll). We have cho- 
sen spectra with a high signal to noise ratio to plot - most of 
the spectra have a significantly lower ratio. 

broad and very shallow, and we have found that minor errors 
in continuum placement can cause a large systematic error in 
the measured centres and widths of these lines. The hydrogen 
Balmer lines were not included in the radial velocity analysis, 
as line blending can cause large and systematic underestima- 
tion of the velocity semiamplitudes (Petrie, Andrews & Scarfe 
119671 Hilditch 1973 Andersen ll975l l. Regions from an exam- 
ple spectrum around the lines used are shown in Fig.|5] 

As there are very few spectral lines strong enough to yield 
reliable radial velocities, we have used four different methods 
to measure these velocities, allowing us to investigate which are 
the most reliable. The techniques are fitting double Gaussian 
functions (hereafter called 2gauss), cross-correlation (onecor), 
two-dimensional cross-correlation (todcor) and spectral disen- 
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Table 3. Additional parameters found during the spectroscopic 
analysis. The full widths at half maximum of the Gaussian 
functions come from the 2gauss analysis and light ratios 
come from the todcor analysis. 



He I line FWsHM of the Gaussian fits ( km s ' ) Light ratio 

(A) Star A Star B from todcor 

4026.189 171.5±9.6 161.9±11.2 0.839+0.138 

4387.928 167.8±10.1 163.8±6.7 0.855+0.106 
4471.681 194.9±16.2 200.8±21.6 0.893+0.059 

4921.929 173.9+9.2 163.8+10.6 0.856+0.082 



tangling. Spectroscopic orbits were calculated for each spectral 
line and for each technique. 

The ONECOR and todcor algorithms require a template spec- 
trum which closely matches the spectral characteristics of the 
components of DW Car. We have used a FEROS spectrum of 
the rotational velocity standard star HR1855 (vOii, spectral 
type BOV), which has a rotational velocity of 1 0-20 km s"' 
(Lyubimkov, Rostopchin & Lambert 120041 Abt, Levato & 
Grosso l2002> . This spectrum has signal to noise ratios of be- 
tween 235 and 350 depending on the echelle order, and was 
broadened by a rotational velocity of 150 km s"' using the 
method of Gray f 1992'\ with linear limb darkening coefficients 
of between 0.39 and 0.30 depending on the wavelength of the 
spectral order. We have measured the radial velocity of this star 
to be H-19.9 + 2.4km s ' from Gaussian fits to several He I 
lines, in agreement with the Wilson-Evans-Batten radial ve- 
locity catalogue (Duflot, Figon & Meyssonnier 1995) value of 
H-17.4kms"'. We have subtracted 19.9kms"' from all the ve- 
locities of DW Car measured using HR 1855 as a template. 

3.1. Gaussian fitting (Igavss) 

The spectra of DWCar were velocity-binned to 5kms ' px ' 
and a double Gaussian function was fitted to four of the five 
good He I lines (listed above) by the method of least squares 
with rejection of points deviating by more than 4o" from an 
initial best fit. The fit to each spectrum was inspected by eye 
and poor fits were rejected. Double weights were assigned to a 
few fits which seemed to be of a particularly high quality. 

Spectroscopic orbits were fitted to the measured radial ve- 
locities for each spectral line and the systemic velocities of 
the two stars were not forced to be equal (see Popper & Hill 
[T99Tt . The orbits were calculated using the method of Lehman- 
Filhes implemented in the sbop^ program (Etzel l2004t . which 
is a modified and expanded version of an earlier code by 
Wolfe, Horak & Storer J1967> . The ephemeris from Paper 1 
was adopted (Table[0 and the orbit was assumed to be circular 
based on initial analyses, the shortness of the orbital period and 
the shape of the light curve. 

The means and standard deviations of the full widths at half 
maximum of the fitted Gaussian functions are given in Table|3] 



^ Spectroscopic Binary Orbit Program, 
jhttp : / /mintaka . sdsu . edu/ f aculty/etzelT] 



3.2. One-dimensional cross-correlation foNEcoRj 

The velocity-binned spectra of DWCar were subjected to a 
cross-correlation analysis (Simkin ri9741 Tonry & Davis [T979l 
using Fast Fourier Transform techniques. We investigated us- 
ing three different types of template spectra to derive radial ve- 
locities from the target spectra: 

- We attempted to use spectrum s2 as a template, which 
was observed at conjunction (phase 0.498), because this 
provides the best match to the target spectra. These at- 
tempts failed because bad pixel values in the target spectra, 
caused by CCD defects, tended to correlate against their 
couterparts in the template spectrum and deform the cross- 
correlation function (CCF). 

- The artificially broadened spectrum of HR 1855 was used 
and found to give good results. This spectrum does not con- 
tain the bad pixels common to the spectra of DW Car We 
therefore used this template to measure radial velocities. 

- Synthetic spectra generated using the uclsyn code (see 
Southworth et al. 2004a for references). The results were 
inaccurate because the synthetic profiles were a poor match 
to the observations, as expected because uclsyn is not in- 
tended for stars this hot. We did not pursue the use of syn- 
thetic spectra further because good results had already been 
obtained using template spectra of HR 1855. 

Radial velocities were accepted or rejected by eye, based 
on the shape of the CCF. The peaks of the CCF were then de- 
termined by least-squares quadratic interpolation and spectro- 
scopic orbits were calculated using sbop as before (Sect. ITU . 
We experimented with using unbroadened template spectra as 
these should be less strongly affected by line blending, but 
found that the results were excessively noisy. This effect has 
been noted before and is due to the template and target spectra 
being quite different to each other. 

3.3. Two-dimensional cross-correlation fxoDcoRj 

The TODCOR algorithm (Zucker & Mazeh ll994> produces a two- 
dimensional CCF by cross-correlating two template spectra 
(one for each star) with each target spectrum, which should re- 
duce the problems caused by blending effects when the stellar 
spectral lines are separated by less than their total broadening. 
As with the onecor analysis, three types of template spectrum 
were investigated and the broadened spectrum of HR 1855 was 
found to be the best. Radial velocities were measured by deter- 
mining the position of the relevant peak of the CCF by fitting 
a surface with the IDL^ procedure min-curve-surf. They were 
accepted or rejected via visual inspection of the CCF, and orbits 
were calculated using sbop as before. 

The TODCOR algorithm can be used to analytically evaluate 
the light ratio between the two component stars (^) on the ba- 
sis of the relative strengths of their absorption lines (Zucker & 
Mazeh 1994 1. The mean and standard deviation of the light ra- 
tio calculated for each spectrum is given in Table |3 for each 
echelle order. These light ratios are meaningful because of the 

^ |http : //www ■ rsinc . com/idTTl 
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Fig. 3. Variation of the quality of the fit for spectral disentan- 
gling the echelle order 55 spectra of DW Car The upper panel 
is a surface plot of fit quality versus Kp, and K^. The lower 
panel is a contour plot of the same data with he six highest val- 
ues of fit quality indicated by filled circles (largest for the best 
fit and smallest for the sixth-best fit). On the lower panel the 
area which a surface was fitted to is indicated by a dashed box 
and the highest point of that surface by a cross. 



negligible difference in spectral type between the two compo- 
nents of DW Car. 



3.4. Spectral disentangling 

The spectral disentangling technique solves simultaneously for 
the contributions of the components to the observed composite 
spectra, and for the Doppler shifts in the component spectra. 
We have applied the method introduced by Simon & Sturm 
J1994> using a version of the original code revised for the 
Linux operating system. It assumes a constant light level, so 
observations taken during eclipse were rejected from the anal- 



ysis; DWCar is constant to within about 5% outside eclipses. 
Dedicated IDL programs were applied to remove cosmic ray 
events and other defects, and for careful normalization of the 
individual orders. 

We have assumed the ephemeris from Paper I and a circular 
orbit, and fitted for the velocity semamplitudes of the spectro- 
scopic orbit for each star, and Kb- Suitable spectral lines are 
present only in orders 55, 51, 50, 45 and 33, and these orders 
were analysed individually. The formal errors in Ka and Kb 
returned by the code are generally around 0.3 km s ', which is 
much smaller than expected given the range in Ka and Kb for 
different orders (Table0] Table|5|l- It is not clear how to calcu- 
late robust uncertainties for spectroscopic orbits calculated in 
disentangling analyses (Hynes & Maxted [l998> . so we will not 
quote uncertainties for individual spectroscopic orbits. 

We have also disentangled the spectra for a grid of fixed 
values of Ka (from 245 to 290km s"') and Kb (260 to 
305 kms ') with a spacing of 1 kms"\ to study how the quality 
of the least-squares fit varies. The results are shown in Fig. |3] 
for the spectra from echelle order 55. There are several local 
maxima but one clear global maximum. The highest point of 
a surface fitted by least squares around the global maximum 
is indicated. This diagram shows that whilst spectral disentan- 
gling can be used to determine spectroscopic orbits for binary 
stars, the best-fitting solution returned by minimisation algo- 
rithms cannot in general be trusted unless confirmed to be the 
global maximum by a grid search or by external information. 

3.5. Corrections from synthetic spectra 

The need to apply corrections to the radial velocities derived 
using the todcor algorithm has been demonstrated in several 
cases, e.g., Torres et al. J1997> and Torres & Ribas ( I2002t . 
These corrections may be needed to remove small radial veloc- 
ity measurement errors due to line blending and/or the shifting 
of spectral features in or out of the analysed wavelength range 
due to orbital motion. The corrections are found by simulating 
observed spectra with known radial velocities and analysing 
them in the same way as in the target spectra. 

A set of simulated spectra were generated for each spectro- 
scopic orbit, using HR 1855 (broadened) as a template, for the 
observed phases and calculated velocities. A light ratio of 0.9 
was used and no additional simulated observational noise was 
added. We have derived these corrections for all four of our ve- 
locity analysis methods because there was a priori no reason to 
assume that the corrections would be negligible for any one of 
them. The derived corrections were applied to the radial veloc- 
ities from the target spectra and new spectroscopic orbits were 
calculated. 

The spectroscopic orbits without velocity corrections are 
given in Table0]and the corrected orbits are given in Table|5j it 
can be seen that in the case of DW Car the velocity corrections 
generally amount to a compensation for the systematic error 
caused by line blending. We have chosen the todcor analysis 
of the He 1 A447 1 line to illustrate these results and the nature of 
the velocity corrections. The uncorrected radial velocities and 
best-fitting spectroscopic orbits are shown in Fig. |3 and the 



Southworth and Clausen: Absolute dimensions of the eclipsing binary DW Carinae 



7 





300 


-1 




1 1 




1 1 1 i_ 

; 




200 


















/ ° 


□\ 












/□ 


□\ 




\ - 
0\ - 


\ 










/o 


o\ \ 


E 


100 


- / 




□ \ 




\ - 

\ 














\ : 















\ : 




r/ 








\^ 


OC 




- \D 






















/ : 


> 














o 


- 1 00 


























o 








o/ 

0/ 




□ / : 














□ / I 




-200 




o/ 








-300 


- 1 




, 1 





- 



0.0 0.2 0.4 0.6 0.8 1.0 

Orbital phase 



Fig. 4. Radial velocities and best-fitting spectroscopic orbit 
measured from analysing the AAAll spectral line with todcor. 
No velocity corrections have been made. Circles represent ve- 
locities for star A and squares for star B. Rejected observations 
are shown using open symbols. They were rejected due to the 
appearance of their CCF, not the value of the derived velocities. 
See text for discussion. 

derived corrections are plotted in Fig. |5] The corrected radial 
velocities are shown in Fig.|6j note that the use of corrections 
allows more velocities to be retained as reliable, improving the 
quality of the resulting spectroscopic orbits. 

The two cross-correlation techniques require the largest 
velocity corrections, which is understandable because CCFs 
are effectively smoothed by the rotational broadening both of 
DWCar and of the template spectrum, which was broadened 
to provide the closest match to the spectra of the components 
of DW Car. Whilst todcor was introduced to minimise the ef- 
fects of spectral line blending (Zucker & Mazeh 1 994) , the 
two-dimensional CCF produced by the todcor algorithm is 
calculated from the normal CCFs between the target and each 
template spectrum and between the two template spectra them- 
selves. As the one-dimensional CCFs suffer from line blending, 
it is not surprising that this effect is still present in the todcor 
CCFs. 

The velocity corrections required by 2gauss were a lot 
smaller than those for onecor or todcor. This is understandable 
because fitting a double Gaussian function explicitly compen- 
sates for spectral line blending. The velocity corrections arise 
mainly due to slight departures in the shape of spectral lines 
from Gaussians. However, the corrections of 2gauss radial ve- 
locities have a larger effect on the resulting systemic velocities 
than those for onecor and todcor. This is again due to slightly 
non-Gaussian spectral line shapes, and is particularly notice- 
able for the /14471 line which has a blended forbidden compo- 
nent on its short- wave length side. 

Velocity corrections to spectroscopic orbits calculated by 
disentangling were negligible, confirming in this case that the 
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Fig. 5. Corrections derived for the radial velocities measured 
from the A447 1 line using todcor. Values for star A are shown 
using filled circles and those for star B using asterisks. These 
corrections are an excellent indicator of the effects of line 
blending. 

technique is an excellent way to measure spectroscopic orbits. 
This suggests that spectral disentangling is not significantly af- 
fected by line blending. This conclusion was also reached by 
Hilditch (|2P05 1. 

3.6. Final Keplerian orbit 

The final velocity semiamplitudes for each radial velocity mea- 
surement technique are given in Table|6lalong with our adopted 
values (velocity semiamplitudes from spectral disentangling 
and systemic velocities from 2gauss). The velocity semiampli- 
tudes are in good agreement with the values derived by Ferrer 
etal. (1985): /STa = 261 + 10kms-i and ^Tb = 278 + lOkms"'. 
Our systemic velocities are in acceptable agreement with the 
-10 + 5 and -16 + 5kms"' found by Ferrer et al. 

We have shown that velocity corrections are needed for 
ONECOR and 2gauss as well as for todcor, and that both cross- 
correlation techniques seem to suffer from line blending even 
when these corrections have been included. We have also dis- 
covered that small differences in how the CCFs are calculated 
(including the size of the spectral range and whether the low- 
est frequency components of the Fourier transforms are sup- 
pressed) can make a significant difference to the uncorrected 
velocity semiamplitudes. These differences normally vanish 
when the velocity corrections are applied, illustrating how im- 
portant the corrections can be. However, strongly blended lines 
continue to give unreliable velocities even when corrections are 
applied. Calculating corrections iteratively, using already cor- 
rected orbits, improves this slightly, but cross-correlation can- 
not reliably deal with spectra which are strongly affected by 
line blending. 

Averaged spectroscopic orbits for each analysis technique 
have been calculated by finding the mean and the standard er- 
ror of the velocity-corrected orbital parameters in Table |5] for 
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Table 4. Spectroscopic orbits calculated for each component of DW Car, given for each of the considered He I lines and for 
each of the four radial velocity measurement techniques. Ka and Kb refer to the velocity semiamplitudes of the stars and Vy^A 
and Vy^B to the systemic velocities (all quantities in kms"'). Systemic velocities cannot be measured as directly using spectral 
disentangling, and no eiTor estimates are available (see text for details). 



Spectral line (A) 


Echelle 
order 




Gaussian fitting 


Cross-correlation 
One-dimensional Two-dimensional 


Spectral 
disentangling 


He I 


55 




259.0 + 2.2 


252.8 ±1.5 


254.5 ± 1.6 


263 


4026.189 




Kb 


277.8 + 1.5 


265.5 ±3.2 


267.3 ±3.3 


277 






Vy.A 


-15.2 ±2.0 


-0.4 ±1.5 


-1.3±1.6 








Vy.B 


-14.5 ±1.3 


-0.7 ±3.1 


-2.3 ±3.2 




He I 


51 


Ka 


261.4 ± 1.3 


259.8 ± 1.8 


257.2 ±1.1 


263 


4387.928 




Kb 


280.6 ±2.0 


280.2 ±2.2 


274.8 ±2.0 


279 






Vy.A 


-3.4 ±0.9 


-1.3 ±1.7 


-1.0±1.0 








Vy.B 


-3.8 ±1.5 


-1.0±2.2 


-0.2 ±1.9 




He I 


50 


Ka 


261.9 ±1.5 


252.2 ±2.0 


255.4 ±1.5 


258 


4471.681 




Kb 


281.5 ±2.2 


274.6 ±3.1 


273.7 ±2.7 


277 






Vy.A 


-30.9 ±1.4 


+9.9 ±2.0 


+0.2 + 1.5 








Vy.B 


-23.7 ±2.0 


+8.3 + 3.1 


+.55 + 2.7 




He I 


45 


Ka 


254.1 ±2.0 


251.6 + 4.0 


251.3 + 3.6 


255 


4921.929 




Kb 


276.3 ±1.1 


269.1 ±3.0 


270.6 + 2.7 


275 






Vy.A 


-11.3±1.8 


-3.4 ±4.0 


-4.4 + 3.4 








Vy.B 


-13.0 ±1.0 


-9.9 ±3.0 


-9.8 + 2.6 




He I 


33 


Ka 








263 


6678.149 




Kb 








281 



Table 5. As Table0]but with velocity corrections applied (see text for details). 



Spectral line (A) 


Echelle 
order 




Gaussian fitting 


Cross- 
One-dimensional 


correlation 

Two-dimensional 


Spectral 
disentangling 


Hel 


55 


Ka 


263.1 +2.1 


262.8 + 1.1 


259.3 + 1.6 


263 


4026.189 




Kb 


282.2+ 1.5 


275.2 + 2.0 


275.0+ 1.8 


277 






Vy.A 


-0.4+ 1.9 


3.6+1.1 


2.0+ 1.5 








Vy.B 


0.4+1.3 


2.5 + 1.9 


0.8+1.6 




Hel 


51 


Ka 


254.6+ 1.0 


251.4+1.2 


252.4+ 1.1 


263 


4387.928 




Kb 


272.1 + 1.7 


266.1 +2.0 


267.6+ 1.8 


279 






Vy.A 


-17.5 + 0.9 


-18.9 + 1.8 


-18.6+1.0 








Vy.B 


-15.5 + 1.5 


-19.2+1.1 


-18.1 + 1.6 




Hel 


50 


Ka 


262.0+1.2 


257.6 + 1.6 


257.8 + 1.4 


258 


4471.681 




Kb 


281.6+1.8 


275.9 + 2.6 


273.9 + 2.2 


277 






Vy.A 


0.0+ 1.1 


5.7 + 2.6 


1.3 + 1.3 








Vy.B 


3.0+ 1.6 


3.0+1.6 


1.4 ±2.0 




Hel 


45 


Ka 


255.0 + 2.0 


252.4 + 2.8 


252.5 ±2.1 


256 


4921.929 




Kb 


276.9 + 1.1 


274.1 + 1.9 


272.2 ± 2.0 


275 






Vy.A 


-6.2+1.5 


-7.4 + 2.6 


-6.9 ± 1.9 








Vy.B 


-8.4+1.1 


-9.2+1.7 


-10.5 ±1.8 




Hel 


33 


Ka 








263 


6678.149 




Kb 








281 



different spectral lines (Table|6j. It can be seen that the velocity 
semiamplitudes calculated using onecor and todcor are lower 
than those from 2gauss and disentangling. The main problem in 
determining these spectroscopic orbits is line blending, which 
will always act to reduce K\ and Kb- Therefore the techniques 
which are least affected by line blending should give higher K^ 
and Kb, implying that disentangling is best, closely followed by 
2gauss, whereas the cross-correlation techniques, onecor and 
TODcoR, are significantly affected. 



We have adopted the averaged spectroscopic orbit (and 
standard error) from disentangling as the best result, with the 
systemic velocities from 2gauss (note that this is not the final 
orbit: see below). The good agreement between this and the av- 
erage orbit for 2gauss shows that the technique of fitting double 
Gaussians also performs very well in situations such as this. 
It is important to note that the results we derived using spec- 
tral disentangling can only be fully trusted because the results 
from other techniques confirm that we have chosen the correct 
peaks in the surfaces of fit quality versus Ka and Kb (Fig.|3}, al- 
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Fig. 6. Radial velocities and best-fitting spectroscopic orbit 
measured from a todcor analysis of the /14471 spectral line 
with with velocity corrections subtracted. Symbols are as in 
Fig.@] The residuals of the fit are shown towards the bottom of 
the plot. Points were rejected on the basis of the apearance of 
their CCF, but taking into account the effects of applying the 
corrections. 



though in every case this peak was a global maximum. We have 
tabulated one set of radial velocities for the interested reader 
(Table0. We chose the 2gauss analysis of the He I A4411 line, 
as 2gauss seems to be more reliable than onecor or todcor, and 
the A441 1 line velocities give a similar result to the adopted fi- 
nal orbit. 

It is also noticeable from Table |6l that the agreement be- 
tween different analyses of the same spectral line is greater than 
that between the same analysis method used on different lines. 
This suggests that the systematic errors of different techniques 
are smaller than the random errors, although it must be remem- 
bered that the systematic errors may remain at full strength 
after the results from several lines have been averaged. The 
systemic velocities from different echelle orders have a poorer 
agreement than expected; any further investigation of this will 
require the use of spectra from a dEB which exhibits a far richer 
set of spectral lines. 

3. 7. Final spectroscopic orbit 

The spectroscopic orbits fitted so far have been purely 
Keplerian, i.e. we have assumed that the velocities of the cen- 
tres of light of the stars are identical to those of the their centres 
of mass. In close binaries, though, the effects of reflection and 
ellipsoidal stellar shape can cause the observed velocities to be 



Table 6. Keplerian spectroscopic orbital parameters for 
DWCar. The mean and standard error are given for each of 
the four radial velocity analysis techniques, and for the final 
results. All quantities are in units of kms"'). 





Ka 


A-B 






V 


2gauss 


258.7 + 1.9 


278.2 + 


2.0 


-7.3 + 3.7 


-6.3 + 4.0 


ONECOR 


256.0 + 2.3 


272.8 + 


2.0 


-4.3 + 4.9 


-5.7 + 4.6 


TODCOR 


255.5 ±1.5 


272.2 + 


1.4 


-5.6 + 4.2 


-6.6 + 4.1 


Spec. dis. 


260.4 ± 1.5 


277.8 + 


0.9 






Adopted 


260.4 + 1.5 


277.8 + 


0.9 


-7.3 + 3.7 


-6.3 + 4.0 



Table 7. Corrected radial velocities from a 2gauss analysis of 
the He I AAAll spectral line. Note that these data are only a 
small part of those used to calculate the final spectroscopic or- 
bits. The columns marked O - C give the observed minus cal- 
culated velocity residuals. 



HJD - 


Weight 


Radial velocities 





-C 


2400000 




star A 


star B 


A 


B 


51144.83060 


1.0 


257.4 


-279.8 


-3.0 


-2.0 


51145.84399 


0.0 


25.0 


13.2 


28.5 


8.8 


51146.86755 


1.0 


-256.0 


284.5 


5.2 


4.5 


51147.86649 


0.0 


-1.4 


48.9 


34.9 


9.4 


51148.80902 


1.0 


247.7 


-278.5 


-12.1 


-1.4 


51149.87170 


0.0 


106.7 


66.2 


55.4 


120.4 


51150.83746 


1.0 


-262.2 


276.5 


0.5 


-5.1 


51152.73657 


1.0 


240.0 


-275.1 


-2.8 


-16.1 


51172.69544 


2.0 


256.2 


-273.6 


-1.2 


1.0 


51174.78088 


2.0 


-252.7 


269.4 


1.3 


-2.8 


51176.80734 


1.0 


246.9 


-245.4 


9.3 


8.0 


51178.84953 


1.0 


-208.1 


226.8 


-1.6 


5.4 


51179.87450 


0.0 


-192.6 


216.3 


-2.2 


12.0 


51180.81349 


1.0 


230.1 


-238.3 


5.8 


1.0 


51181.87490 


1.0 


193.7 


-219.9 


-2.4 


-10.8 


51182.86220 


0.0 


-184.0 


209.9 


-2.2 


14.8 


51184.86864 


0.0 


167.6 


-193.0 


1.1 


-15.6 


51187.87364 


1.0 


-220.4 


233.8 


-4.4 


2.1 


51189.86682 


1.0 


215.1 


-223.6 


-0.5 


6.2 


51190.84867 


0.0 


-207.3 


118.5 


-44.1 


-56.6 


51191.86279 


1.0 


-229.4 


238.5 


-9.3 


2.5 


51192.80846 


1.0 


199.0 


-201.8 


8.4 


1.4 


51193.87953 


1.0 


232.0 


-244.4 


-2.1 


5.2 


51197.85496 


1.0 


230.3 


-254.2 


0.7 


-9.4 


51199.87004 


2.0 


-247.1 


254.8 


-3.5 


-6.4 


51209.72395 


1.0 


164.3 


-196.2 


-1.9 


-19.1 


51382.51044 


2.0 


263.2 


-270.3 


1.9 


8.6 


51386.50761 


2.0 


261.8 


-273.0 


0.5 


5.9 


51563.86655 


2.0 


-221.5 


251.4 


6.1 


7.4 



Table 8. Final spectroscopic orbit after correction for non- 
Keplerian effects. / and a are the orbital inclination and semi- 
major axis, respectively. 



Star A 



StarB 



K (Keplerian) ( km s" ) 
K (corrected) (kms"') 
Mass ratio 
a sin i ( Rq) 
M sin''/ (Mo) 



260.4+ 1.5 277.8 + 0.9 
261.7 + 1.6 279.3 + 1.2 
0.9370 + 0.0070 
14.192 + 0.052 
11.25+0.12 10.54 + 0.14 
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significantly different to the centre-of-mass velocities. These 
effects are accurately modelled by the Wilson-Devinney code 
(see Section which can be used to correct individual ra- 
dial velocities using a given system geometry. However, our 
adopted velocity semiamplitudes were calculated using spec- 
tral disentangling. With this technique the spectroscopic orbit 
is calculated directly without the intermediate step of calculat- 
ing radial velocities. We therefore cannot apply radial velocity 
corrections and so must calculate corrections which can be ap- 
plied to the velocity semiamplitudes themselves. 

We have calculated spectroscopic orbits for the 2gauss ra- 
dial velocities, using the Wilson-Devinney code both with and 
without the inclusion of non-Keplerian effects. This has been 
done for echelle orders 55, 51, 50 and 45. The mean value 
and standard error of the resulting corrections to the veloc- 
ity semiamplitudes (in the sense Ktrua - ^^Kepierian - ^K) are 
AKa = -Hl.3 + 0.1kms-i and A/Tb = -h1.5 + 0.2kms-i. This 
indicates that in this case fitting Keplerian orbits caused the 
velocity amplitudes (and so the masses) to be systematically 
underestimated by a small but significant amount when fitting 
a Keplerian orbit. This is the expected result, because the light 
centres for each star are shifted towards the companion due to 
their mutual irradiation, causing the measured radial velocities 
to be underestimates of the motion of the centre of mass. 

The final spectroscopic orbit, with the non-Keplerian cor- 
rection applied, is given in Table |S1 The resulting minimum 
masses are M^sin^ i - 11.25 + 0.12 and Mb sin"' / = 
10.54 + 0.14 Mq where / is the orbital inclination. The correc- 
tions to the systemic velocities are negligible and have not been 
applied. 

3.8. Spectroscopic light ratio 

The light ratio of the two stars is poorly defined by the light 
curves of the system (see below) so we have derived a spec- 
troscopic light ratio to help constrain the light curve solutions. 
The TODCOR algorithm was used to analytically determine the 
light ratio for each spectral line analysed (see Sect. 13.31 and 
Table |3j, yielding an average value and standard deviation of 
^ = 0.873 + 0.042. The wavelength dependence of the light 

ratio is negligible between the Balmer jump and 5500 A (see 
Section Rm . This is a reliable equivalent width ratio measure- 
ment because the components of DWCar have very similar 
spectral characteristics, and the same template spectrum was 
used for both stars. 

The equivalent widths of the He 1 A4026 and A441 1 lines 
for both stars were also measured directly using several spec- 
tra with large velocity separations and high signal to noise ra- 
tios. The component lines are still slightly blended so we have 
used the deblending feature of the iraf^ task splot to measure 
the line strengths. The mean values and standard deviations are 
0.87 + 0.03 and 0.89 + 0.03 for the ^4026 and ^4471 lines re- 
spectively, in good agreement with the todcor results. 

"* iRAF is distributed by the National Optical Astronomical 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc. under contract with the National 
Science Foundation. 



The He I lines are intrinsically slightly stronger in star B 
(Feff = 26 500 K) than in star A (reff = 27 900 K), so the equiv- 
alent width ratios are slightly different to light ratios. A correc- 
tion was derived using intrinsic equivalent width ratios taken 
from ucLSYN and was found to be 0.97 + 0.01. The final light 
ratio, using the equivalent width measurements, is 0.84 + 0.04 
for the /14026 line, which is close to the central wavelength of 
the Stromgren v passband, 41 10 A (Stromgren ll963> . 

3.9. Spectroscopic helium abundance 

At the suggestion of the referee we have attempted to derive 
surface helium abundances for the components of DWCar 
The high rotational velocities mean that this is difficult, but 
spectroscopic abundance analyses have been successfully per- 
formed on similar objects by Pavlovski & Hensberge J2005> 
and Pavlovski et al. (200^. Our data contain insufficient infor- 
mation to significantly constrain the metal abundances of the 
components of DW Car. 

Model atmospheres were calculated using the atlas9 code 
(Kurucz [1979> . Non-LTE atomic level populations were ob- 
tained for hydrogen and helium using the detah. code and syn- 
thetic spectra were constructed using the surface code (Butler 
& Giddings'1985). We adopted the effective temperatures and 
light ratio for DW Car found below and calculated the residuals 
of the best fit for several different helium abundances. 

Formally, the best fits are found for a helium abundance 
equal to or slightly lower than the solar value (Fig.0, but with 
a large uncertainty. Varying the metal abundance within reason- 
able limits has no significant effect. However, whilst the fit of 
the lines at 4388 A and 447 1 A are excellent, the 4026 A fine is 
too shallow in the synthetic spectrum whilst the 4922 A line is 
too deep. The investigation of these effects is beyond the scope 
of this work. We conclude that the spectra are consistent with a 
solar or subsolar helium abundance but not with an abundance 
signficantly greater than solar 

The rotational velocities were also optimised when fitting 
the theoretical spectra to the disentangled spectra, and the best- 
fitting values are Va sin / = 182 + 3 km s"' and Vb sin / = 177 ± 
3kms"^ 

4. Light curve analysis 

The light curves of DWCar contain 518 photoelectric observa- 
tions in each of the uvby passbands (Paper 1). The two eclipses 
are relatively long in phase units and are both about 0.7 mag 
deep, indicating that the system is composed of two similar 
and close stars. We have used the 2003 version (Van Hamme 
& Wilson 2003 1 of the Wilson-Devinney code (Wilson & 
Devinney 1971 Wilson [D7^ll993> . hereafter called wd2003, 
to model the light curves, as Roche geometry is needed to 
adequately represent the tidal deformation of the two stars. 
wd2003 optionally uses the predictions of the atlas9 model at- 
mospheres (Kurucz 1993i) to connect the flux ratios of the two 
stars in different passbands. 

Initial investigations showed that the orbital eccentricity is 
negligible, so we assumed that the orbit is circular. The rota- 
tional velocities of the stars were fixed at the synchronous val- 
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Fig. 7. Plot of the normalised disentangled profiles of four helium lines for each star compared to theoretical spectra calculated 
for helium abundance Y = 0.28 and the effective temperatures and light ratio found in our analysis. The left-hand panels are for 
star A and the right-hand panels for star B. The observed spectra are plotted using points and the theoretical spectra using solid 
lines. 



ues as these are consistent with the spectral line widths. The 
bolometric albedo coefficients and gravity brightening expo- 
nents were fixed at 1.0, which is the expected value for stars 
with radiative envelopes (Claret I1998> . The mass ratio was 
fixed at the spectroscopic result of 0.937 and we chose to model 
reflection using the faster and simpler of the two alternatives 
implemented in wd2003. We have used grid sizes of = 20 
for each star; solutions with finer grids gave negligibly differ- 
ent results. 

We expected that there would be a small amount of contam- 
inating 'third' light in the light curves of DW Car, as its Be star 
nature indicates that there is substantial circumbinary material. 
Third light was therefore included as a fitting parameter in the 
light curve analysis, but the optimised values were all consis- 
tent with, and marginally lower than, zero. We therefore have 
assumed that third light is zero. A smaU amount of third light 
(up to about 5%, which is an order of magnitude greater than 
the best-fitting values) will have a negligible effect on the light 
curve solutions. 

Limb darkening coefficients were initially set to values 
found by bilinear interpolation in effective temperature and sur- 
face gravity from the coefficients tabulated by Van Hamme 
J1993> . but in our final solutions we have always optimised 
these values to avoid the theoretical dependence caused by us- 
ing coefficients calculated from stellar atmosphere models. The 



choice of linear, logarithmic or square-root law has a negligi- 
ble effect on the solutions so we have chosen to use the latter, 
as this provides the best fit to model atmosphere predictions 
for hot stars (Van Hamme I1993> . Fixing the limb darkening 
coefficients at either the Van Hamme (I 1993> or Claret ( I2000> 
theoretical values causes a negligible change in the final solu- 
tion - of the order of 0.1% in the stellar radii. The optimised 
values are generally in good agreement with the theoretically 
predicted coefficients. 

The radii found by fitting the light curves with wd2003 are 
strongly correlated with each other and with the light ratio, re- 
sulting in one of many different light curve solutions (i.e., local 
minima in parameter space) being found depending on the ini- 
tial fitting estimates. The residuals of the fit are approximately 
constant for ratios of the radii between 0.8 and 1.2. We have 
therefore used the light ratio of 0.84 + 0.04 found in Sect.lT^to 
constrain the solution of the v band light curve. This was done 
by fixing the potential of star A at values which gave the de- 
sired light ratio and plus or minus its uncertainty. These values 
for potential were then fixed for the uby light curve solutions. 
The potential of star B, the orbital inclination, the phase shift 
between primary midminimum and phase zero, and the light 
contributions for the two stars were optimised for each light 
curve. The linear terms in the square-root limb darkening law 
were also optimised. The effective temperatures of the two stars 
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were fixed at nominal values (27 500 and 26 750 K), but these 
do not affect the final solution because we are fitting for the 
light contributions of both stars directly. Fig. |8l shows the light 
curve fits and residuals. 



Table |5] contains the values of the optimised parameters, as 
well as the standard errors calculated by wd2003, for each light 
curve and for each of the three potentials of star A (4. 101, 4. 132 
and 4. 157) which are required to reproduce the light ratios 0.80, 
0.84 and 0.88 for the v light curve. We have adopted the results 
for a star A potential of 4.132 as the final values. The uncertain- 
ties in these light curve parameter values have been calculated 
by adding in quadrature an uncertainty from the variation for 
different potentials to the mean standard error for each param- 
eter. Additional uncertainties due to the treatment of third light 
and limb darkening are negligible. 



4.1. Stromgren indices for each star 

Separate Stromgren photometric indices can be derived for 
each star from the indices of the system and the light ratios in 
the uvby passbands. This requires a solution of the light curve 
with a consistent geometry, so we have derived the flux ratios 
and corresponding Stromgren indices (and uncertainties) by fit- 
ting the light curves with the potential of star A fixed at the 
values found above. For each fixed potential, we have calcu- 
lated the Stromgren indices for each star from the indices of 
the combined system and the light ratios in the four light curves 
(Table|9j. The uncertainties in TableOhave been calculated by 
adding in quadrature all the relevant uncertainties, including 
those resulting from the range of possible values of the poten- 
tial of star A. They are somewhat lower than they would be if 
we had simply used the final light ratio errors for each light 
curve, because the light ratios are correlated with each other 
(i.e., a change in the potential of star A causes absolute changes 



Southworth and Clausen: Absolute dimensions of the eclipsing binary DW Carinae 



13 



Table 9. wd2003 light curve solutions for DW Car. A light ratio was used to constrain the solution for the v light curve, and the 
resulting potential of star A was fixed for the other light curves (see text for further details). The quoted errors (except for the 
adopted values in the final column) are the standard errors calculated by wd2003. 



Parameter 


Potential 


V 


y 


b 


u 


Adopted 




of star A 


+ 


+ 


± 


± 


values 


Light ratio 


4. 


101 


0.800 


0.796 ±0.012 


0.797 ± 0.009 


0.769 ±0.034 






4. 


132 


U.o4U 


U.o4U ± U.U14 


U.o4U ± U.UiU 


U.ol J ± U.Uo4 






/I 


1 ^7 


0.880 


870 + on 


SSI +0010 


SSS + 0^7 




Potential of 


4. 


101 


4.225 ± 0.007 


4.236 ± 0.009 


4.228 ± 0.007 


4.248 + 0.011 




star B 


4. 


132 


4. 1 VZ ± U.UU / 


4.ZUU ± U.UuV 


A 1 -u n nn7 
4. IV J + u.uu / 


/I 111 -u n ni n 
4.Z1 i ± u.uiu 






4. 




4 1 ■S7 + 007 


A \fn + OOS 


4 1 "SS + 007 


A 177+0 OOQ 




Fractional radius 


4. 


,101 


0.3237 


0.3236 


0.3237 


0.3237 


0.3203 


of star A 


4. 


,132 


U. j/Uj 


U.jZUj 


v.ilKji 




+ U.UUjU 




4. 


1 -^7 


0.3176 


0.3175 


0.3175 


0.3175 




Fractional radius 


4. 


,101 


0.2993 ± 0.0007 


0.2982 ± 0.0009 


0.2990 ±0.0007 


0.2970 ±0.0011 


0.3019 


of star B 


4. 


,132 


U. JUZO ± u.uuu / 


n 1 s -1- n nnriQ 
U.JUlo ± U.UUUV 


u.JUzj ± u.uuu/ 


U.JUU/ ±U.UUiU 


+ U.UUj / 




A 

4. 


1 "^7 


3063 + 0007 


^052 + 0009 


3061 + 0007 


3048 + 00 1 




Orbital inclination 


4. 


,101 


85.76 ±0.07 


86.00 ± 0.09 


85.96 ± 0.07 


86.13 ±0.08 


85.72 


(degrees) 


4. 


,132 


oj.ol ± U.U/ 


Qj. 11+ U.U4 


OJ. /j + U.U4 


oj. /4 ± U.Uo 


+ U.Zo 




A 

4. 


1 '^7 




7^ + OS 


SS 70 + 07 


Q-^ c:i + OQ 




Phase shift of primary 


4. 


,101 


0.06 ± 0.04 


-0.02 + 0.03 


-0.02 ± 0.05 


0.03+0.05 




midminimum (10 ^) 


4. 


,132 


U.Uo ± U.Uj 




-U.Ui + U.Uj 


U.Ui + U.Uj 






A 

4. 


1 "^7 




_n 04 + 04 


00 + OS 


04 + 04 




Linear LD coefficient 


4. 


101 


-0.124 ±0.031 


-0.132 ±0.037 


-0.156 ±0.029 


-0.056 ± 0.038 




of star A (using the 


4. 


,132 


r\ no /I _i_ r\ f\i i 
-U.Uo4 ± U.Ujl 


-U.Uy 1 ± U.Ujo 


-U. 1 iJ ± U.UjU 


n mi _i_ n m o 
-U.UZj ± U.Ujo 




square-root law) 


4. 


157 


-0.125 ±0.036 


-0.1 14 ±0.043 


-0.137 ±0.036 


-0.039 ± 0.045 




Linear LD coefficient 


4. 


101 


-0.093 ± 0.042 


-0.087 ±0.050 


-0.1 17 ±0.041 


-0.069 ±0.058 




of star B (using the 


4. 


132 


-0.095 ± 0.042 


-0.068 ±0.052 


-0.109 ±0.042 


-0.017 ±0.057 




square-root law) 


4. 


157 


-0.145 ±0.043 


-0.1 14 ±0.052 


-0.149 ±0.042 


-0.065 ±0.058 




Size of residuals 


4. 


101 


5.21 


6.20 


4.97 


7.15 




CTrms (mmag) 


4. 


,132 


5.16 


6.23 


5.00 


7.04 






4. 


,157 


5.39 


6.23 


5.02 


7.04 





in the light ratios significantly larger than they change relative 
to each other). 

The uvhy light curves have also been solved simultane- 
ously to derive an effective temperature difference for the two 
component stars. In this case, the light ratios in different pass- 
bands must be connected, which can be achieved either using 
Planck law or by using the predictions of theoretical model at- 
mospheres. The Planck law is not a good option for DWCar 
because the Balmer jump causes the stars' radiative proper- 
ties to be quite different from black bodies. wd2003 option- 
ally uses predictions from Kurucz atlas9 model atmospheres, 
which should be reliable in this case as the effective tem- 
peratures of the two stars are very similar. Using this op- 
tion, we have obtained an effective temperature difference of 
TeffA - TeffB = 1380 + 420 K for the stars in DWCar, where 
the error has been found by adding in quadrature the uncer- 
tainty caused by varying the potential of star A as before, and 
the formal error of the temperature difference calculated from 
the covariance matrix by wd2003. Formal errors can be un- 
derestimates of the real uncertainties of light curve parameters 
(Popper 1 1984t Maceroni & Rucinski [T997l Southworth et al. 



Table 10. Light ratios of DW Car for a consistent geometry, 
with the resulting Stromgren indices. Note that the quoted flux 
ratio errors were not used in calculating the errors of the indices 
of each star. 







y 


h 




u 


Flux ratios 




0.840 


0.840 


0.840 


0.813 




± 


0.046 


0.044 


0.040 


0.056 








V 


h-y 


m\ 


C\ 


System 




9.675 


0.067 


0.046 


-0.017 




± 


0.005 


0.005 


0.008 


0.008 


Star A 




10.307 


0.067 


0.046 


-0.034 




± 


0.026 


0.011 


0.016 


0.023 


StarB 




10.563 


0.067 


0.046 


0.004 




± 


0.032 


0.014 


0.020 


0.030 



5. Absolute dimensions 

The masses, radii and surface gravities of the two components 
of DWCar have been calculated from the results of the spec- 
troscopic and photometric analyses. Using careful error propa- 
gation, we find that the masses and radii are known to accura- 
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Table 11. The physical parameters of the component stars of 
DW Car derived from the results of the spectroscopic and pho- 
tometric analyses. Vgq denotes the observed equatorial rota- 
tional velocities of the stars. 

t Calculated assuming Lg = 3.826x10^*' W and Mboi© = 4.75 
(Zombeck 1990). 



Parameter 


Star A 


StarB 


Mass (Mq) 


11.34±0.12 


10.63 + 0.14 


Radius (Rq) 


4.558 + 0.045 


4.297 + 0.055 


logfe/cm s"3) 


4.175 ±0.008 


4.198 + 0.011 


Veq(kms-') 


183 + 3 


178 + 3 


Vsyt,ch (kms"') 


173.6 + 1.7 


163.7 + 2.1 


logTeff (K) 


4.446 + 0.016 


4.423 + 0.016 


log(L/Lo)t 


4.055 + 0.063 


3.915 + 0.067 




-5.39 + 0.16 


-5.04 + 0.17 



cies of about 1% (Table fTTl . DWCar therefore joins the (very 
short) list of well-studied high-mass unevolved dEBs, along- 
side QX Car (Andersen et al. ll983> and V578 Mon (Hensberge, 
Pavlovski & Verschueren l2000l l. 

We have derived the effective temperatures of the two stars 
from the Stromgren indices found in Section 14.11 using the 
calibration of Moon & Dworetsky ( 1985) contained in the 
code uvBYBETA (Moou 1985 ) . Using the individual Stromgren 
indices for each star, we find TeffA = 27 800 + 1100 and 
TeffB = 26 500 ± 1300K, where the error does not include any 
uncertainty in the calibration and its main contribution comes 
from the uncertainty in the mi indices. However, the equiva- 
lent effective temperature of the system is much better defined: 
T^eff A+B = 27 230 + 390 K, and the temperature difference is 
also well defined from the light curve analysis: Teff a - T^effB = 
1380+420K (formal error). Using these quantities, and includ- 
ing an estimated uncertainty in the photometric calibration, we 
find temperatures of 27 900 + 1000 and 26 500 + 1000 for the 
two stars. The calibrations of Davis & Shobbrook ( 1977) and 
Napiwotzki, Schonberner & Wenske ( 1993 > give temperatures 
in good agreement with these values. 

The equatorial rotational velocities quoted for each star 
were found by comparing the disentangled spectra of the stars 
to theoretical spectra in order to determine the surface helium 
abundance (Section ll!9l . They are slightly greater than the syn- 
chronous value despite the short orbital period of DW Car. This 
indicates that it is a young system in which the rotational and 
orbital velocities have not yet been equalised by tidal effects. 

The reddening of the system has been found using the 
UVBYBETA codc, which implements the {b - y)o vs. cq rela- 
tion by Crawford ( il978l ). The reddening for the system in- 
dices and for the individual stars is = 0.180 + 0.007, 
Eb^y ^ = 0.177 + 0.054 and Et-y g = 0.182 ± 0.064, where 
the quoted errors come only from varying the indices within 
their uncertainties. The main contribution to the reddening un- 
certainties for the individual stars comes from the uncertainty 
in their {b - y) colours. We will adopt Eh-y = 0.18 + 0.02 for 
DW Car, where the quoted error includes a contribution from 
the estimated uncertainty in the calibration. 



5.1. The distance to DWCar 

Now we have the absolute dimensions and effective temper- 
atures of the two stars we can calculate the distance to the 
system. Distances are in general more accurate and precise 
at infrared wavelengths (Southworth et al. "2005a'), but we are 
lacking reliable apparent magnitudes in the JHKL passbands. 
These would be available from 2MASS, but are too faint be- 
cause they were taken during secondary eclipse. JHKL pho- 
tometry is available from Tapia et al. ( 1988) but seems to be 
too bright compared to the optical colours of DWCar. The 
emission-line nature of DWCar suggests that there may be 
a circumbinary disc. This may affect the infrared colours of 
the system, resulting in brighter apparent magnitudes. We are 
therefore unable to use infrared magnitudes of DW Car to mea- 
sure its distance with any confidence. 

We have used the V band magnitude and flux ratio of 
DW Car (Table I10> to determine the distance to each star and 
to the system using several sources of bolometric corrections. 
We used the jktabsdim code'', which calculates distances using 
several different sources of bolometric corrections (Southworth 
et al. 2005a). The available sets of bolometric corrections are 
calculated using, in general, different radiative parameters for 
the Sun (Bessefl, Castefli & Plez 1998 ). It is important to adopt 
the appropriate values when using bolometric corrections, and 
improper results can be obtained if this is not done (see discus- 
sion in Southworth et al. l2005cl . 

Adopting the bolometric correction formula of Flower 
(I1996i . we find that the individual distances of the two stars 
are in acceptable agreement (2877 and 2742 pc), and the over- 
all distance measurement to the system is 28 10 + 160 pc, where 
the main contributor to the uncertainty comes from the redden- 
ing. If we use the theoretically-derived V band bolometric cor- 
rections of Bessefl et al. (1998) or Girardi et al., (2002) we find 
distances to DW Car of 2710+ 140 pc and 2720+ 140 pc, respec- 
tively. The agreement between the different distance measure- 
ments is acceptable given the uncertainties, but not as good as 
we were expecting. The components of DWCar are hot stars, 
and transforming their visual absolute magnitudes into bolo- 
metric ones requires corrections which are large and relatively 
uncertain. An improved distance estimate could be found by 
studying the ultraviolet to optical spectral energy distribution 
of the system (see Guinan et al. ll998l and Fitzpatrick & Massa 
.1999J . 

Using the RI photometry of DW Car from Carraro & Patat 
(1200 U and bolometric corrections from Girardi et al. (I2002> 
gives distances of 2890 ± 130pc and 2710 + 1 lOpc, for R and 
/ respectively, where the decreasing importance of the uncer- 
tainty in reddening towards longer wavelengths can be seen. 
However, we will adopt 2810 ± 160pc as our final distance, 
found from the V band photometry and the Flower bolometric 
correction formula, because this is an empirical measurement. 
The corresponding distance modulus is 12.24 ± 0.13 mag. 



^ This can be obtained from 
http : //www, astro . keele . ac . uk:~ jkt/co^es .htmll 
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6. Comparison with theoretical models 

Graphical comparisons between the predictions of theoretical 
stellar models and the physical properties of a dEB provide a 
clear and straightforward interpretation of the agreement or dis- 
agreement between theory and observation. The most impor- 
tant plot is of mass versus radius, as these are the two quantities 
which have been directly measured for the dEB (i.e., without 
using any calibrations). A second plot can be made because we 
have measured both the radii and effective temperatures of two 
stars with known masses. Chosing the ordinate quantity to be 
mass is again the best approach because the masses of the stars 
have been directly measured and because model predictions are 
calculated for specific (initial) masses. The abscissal quantity 
could be either effective temperature or luminosity (or, equiv- 
alently, absolute bolometric magnitude); we prefer the former 
because this quantity has been measured more directly than lu- 
minosity. 

In dEBs such as DW Car, which are composed of two simi- 
lar stars, the effective temperatures can be very similar and the 
errorbars often overlap. In this case, the temperature difference 
suggested by the individual temperatures (Table [TTt is 1400 + 
1420 K, whereas results from the light curve analysis give a 
much more precise temperature difference of 1380 + 420 K. 
This can be represented on a mass-temperature plot by giving 
full-size errorbars to one star (here star A) and only relative 
errorbars to the other star (B). 

We have chosen to compare the properties of DW Car with 
theoretical predictions from the Claret evolutionary models 
(Claret 1 995. ,1997. Claret & Gimenez il995 j which are avail- 
able for several fractional helium abundances, Y, for each frac- 
tional metal abundance, Z. These models use the opacities of 
Iglesias, Rogers & Wilson (.1992.) and incorporate moderate 
convective core overshooting (with aov = 0.20) and mass loss, 
and are given for scaled-solar metal abundances. 

Isochrones for the best-fitting ages are plotted in Fig.|9]for 
metal abundances of Z = 0.004, 0.01 and 0.02. A subsolar 
metal abundance, Z = 0.01, provides an excellent match to 
the masses and radii of DW Car, but predictions for the other 
chemical compositions are also in reasonable agreement here. 
However, the standard helium abundances in the Claret models 
(dashed lines) give rise to predicted temperatures slightly too 
cool for Z - 0.01 and 0.02, although good agreement is found 
for Z = 0.004. A slightly enhanced helium abundance gives 
a better fit. All models predict shallower mass-T^eff slope for 
DWCar, although they agree within the uncertainties. 

Other sets of model predictions (Schaller et al. 1992 Pols 
et al. 1199 8 ) were investigated and similar results were found. 
'Classical' models (without convective core overshooting; Pols 
et al. I1998> . gave predictions which were negligibly differ- 
ent from their equivalent overshooting models, as expected for 
such young stars. 

Whilst we found the surface helium abundance of DW Car 
to be approximately solar (Section f3.9> . the model comparison 
suggests an overall abundance slightly greater than solar Note 
that these two statements are not mutually inconsistent because 
they refer to different parts of the stars. However, neither abun- 
dance measurement is very accurate. We conclude that the best 



match between stellar models and the properties of DWCar is 
found for a metal abundance of Z = 0.01, a helium abundance 
of F = 0.26 (roughly solar), and an age of 6Myr. However, 
these three quantities are all quite uncertain and have a depen- 
dence on theoretical models. The age and chemical composi- 
tion we infer for DW Car are relatively imprecise because its 
two component stars have similar properties. Much more pre- 
cise conclusions are possible for dEBs with smaller mass ratios 
(e.g., Southworth et al. l2004bl . 



7. DWCar and Collinder228 

DW Car is a photometric member of Cr 228 and has a sky po- 
sition close to the cluster centre. However, the systemic veloc- 
ities of Cr228 quoted by Rastorguev et al. ( 1999 ) and Levato 
et al. ( I1990> . -12kms ' and -13.5kms ' respectively, dis- 
agree slightly with our measured systemic velocity for DW Car 
(Table |6}. However, as noted in Sect. 13.71 the systemic ve- 
locities from the different echelle orders are not in full agree- 
ment considering their formal errors. Ferrer et al. (I1985> found 
systemic velocities for the components of DWCar in good 
agreement with the cluster velocity. As our measured distance 
and reddening are also appropriate for cluster membership, it 
seems most probable that DWCar does belong to Cr228. We 
can therefore derive the important astrophysical properties for 
Cr 228 from those for DW Car. 

The distance we derive, (m - M)o - 12.24 + 0.13 mag, is 
in good agreement with most of the measurements in the lit- 
erature, including that of The et al. (I1980> for a normal red- 
dening law, but not including those of Carraro & Patat ( 1200 U 
and Tapia et al. (.1988 ). The disagreement with Tapia et al. re- 
sults from the very large interstellar extinction value they find 
using infrared photometry. Our distance is the most accurate of 
the measurements, demonstrating that dEBs are excellent dis- 
tance indicators. Our reddening value for DWCar, equivalent 
to Eb-v = 0.25 + 0.03 mag (Crawford '1978\ is in reasonable 
agreement with most literature values. 

The age of Cr 228 is less well known. Determinations on 
the basis of photometry alone are tricky due to the sparse na- 
ture of the cluster, difficulties with measuring reddening, and 
because UBV photometry is a poor temperature indicator for 
hot stars (e.g., Massey & Johnson [1993}. Massey et al. (1200 U 
used classification spectroscopy to measure ages for many stars 
in Cr 228 and found an age of about 2 Myr, in good agreement 
with our estimate of about 5 Myr from comparison with the 
Claret evolutionary models. 

The chemical composition of Cr228 has not previously 
been studied. We find a metal abundance of Z ^ 0.01 and 
an approximately solar helium abundance from the properties 
of DW Car A spectroscopic abundance study of some of the 
more slowly-rotating members of Cr 228 would be very useful 
in specifying the chemical composition of the cluster and of 
DWCar, which would provide an additional constraint when 
comparing model predictions to the properties of this system. 
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Fig, 9. Comparison between the properties of DWCar and the predictions of the Claret stellar evolutionary models. Panels (a) 
and (b) are a simultaneous match for metal abundance Z = 0.004, panels (c) and (d) likewise for Z = 0.01 and panels (e) and (f) 
for Z = 0.02. The model helium abundances and interpolated ages are given in the legend of each left-hand panel. The errorbars 
plotted for the effective temperature of star B are relative to the effective temperature of star A and do not represent the overall 
uncertainty. 
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8. Summary 

DW Carinae is a young, high-mass detached eclipsing binary 
system with a short orbital period. It is particularly interesting 
because it is unevolved, shows Be star emission-line character- 
istics, and is a member of the young open cluster CoUinder 228 
in the Carina complex. 

The spectra of DW Car contain strong hydrogen lines, with 
some superimposed Be-type emission from circumbinary ma- 
terial, but there are very few other noticable absorption lines. 
As a result of this, we were only able to measure radial ve- 
locities from four or five helium lines. The high rotational ve- 
locities of the stars mean that other lines are not detectable, or 
that small changes in continuum placement can cause a large 
change in the derived radial velocities. Due to the high rota- 
tional velocities, we measured radial velocities using several 
ways to investigate which are the most reUable in the presence 
of strong line blending. 

Individual radial velocities were obtained by fitting double 
Gaussians, by cross-correlation, and by the two-dimensional 
cross-correlation algorithm todcor (Zucker & Mazeh 1994 1 us- 
ing the spectrum of a standard star as a template. Spectroscopic 
orbits were also obtained directly (by least-squares fitting) us- 
ing spectral disentangling (Simon & Sturm .1994 ). although 
this procedure is not straightforward due to a profusion of lo- 
cal maxima in the surface of the quality of the fit in parameter 
space. 

We have found that the radial velocities derived from 
Gaussian fitting and from cross-correlation require substan- 
tial corrections for the effects of line blending. These cor- 
rections were derived from simulated composite spectra, con- 
structed using a broadened observed template spectrum, for ra- 
dial velocity offsets appropriate for each observed spectrum. 
The corrections were the best-behaved for Gaussian fitting, de- 
spite the obvious shortcoming of this technique that helium 
lines do not have strictly Gaussian shapes. However, the cor- 
rections for cross-correlation do not seem to have been com- 
pletely successful, in that the resulting spectroscopic orbits 
have systematically lower ampitudes that those derived from 
the other two methods. The line blending problem is more se- 
vere for cross-correlation, because both the object and the tem- 
plate spectra are broadened (using an unbroadened template 
spectrum causes excessive random errors in the resulting ra- 
dial velocities). Gaussian fitting, on the other hand, explicitly 
accounts for the presence of blending lines from the other star, 
so gives results which are reliable. The cross-correlation results 
could probably be improved by fitting single-lined CCFs to the 
double-lined CCFs (see Hill ll993j . The broadening-function 
approach advocated by Rucinski (I2002> is also a useful alter- 
native. The poor performance of todcor in dealing with line 
blending arises because it presents no significant advantage 
over traditional cross-correlation when the component stars 
have very similar spectral characteristics (P. F. L. Maxted, pri- 
vate communication). 

The orbits measured from spectral disentangling were the 
most internally consistent, and did not require correcting for 
blending. Disentangling should be innately more appropriate 
to this situation because it assumes nothing about the spectral 



characteristics of the stars, and uses (almost) every observed 
spectrum in a simultaneous least-squares fit. However, it suf- 
fers from the presence of many local maxima in the surface of 
quality of the fit to the observations. The best approach for mea- 
suring spectroscopic orbits in the presence of strong line blend- 
ing is therefore to use spectral disentangling, checked and cor- 
roborated by the results of fitting double Gaussian functions. 
Conducting such analyses for many lines (helped if the spectra 
have the wide wavelength coverage achievable using echelle 
spectrographs) is very useful as it allows the overall errors to 
be both reduced and robustly estimated. 

The light curves of DWCar were modelled using the 
Wilson-Devinney code. The ratio of the radii of the stars is 
intrinsically poorly determined by the light curves due to cor- 
relations between various photometric parameters, as is often 
found for systems with deep but not total eclipses. This degen- 
eracy was broken by constraining the light curve solution using 
a spectroscopically-derived light ratio. We also fitted directly 
for the limb darkening coefficients of the stars in order to avoid 
a dependence on theoretically-derived coefficients. The adop- 
tion of different limb darkening laws (linear, square-root and 
logarithmic) makes a negligible difference to the result, and the 
fitted coefficients are also in good agreement with the predic- 
tions from model atmosphere analyses. 

The masses and radii of the components of DW Car, cal- 
culated from the results of the spectroscopic and photometric 
analyses, are known to accuracies of about 1%, so this sys- 
tem joins the short list of well-studied dEBs containing com- 
ponents with masses above 10 Mq. The effective temperatures, 
and interstellar extinction were derived using photometric cali- 
brations and the individual Stromgren indices, which were cal- 
culated for the two stars from the combined indices and the 
light ratios found during the hght curve analysis. 

The properties of DW Car can be used to infer the proper- 
ties of its parent cluster, Cr228. The distance modulus of the 
dEB, calculated using the empirical bolometric corrections of 
Flower J1996> . is 12.24 + 0.13 mag, which is the most accurate 
measurement of Cr228's distance yet. A comparison between 
the properties of DW Car and predictions from the theoretical 
evolutionary models of Claret ( 1995 1 was also undertaken. The 
mass-radius and mass-Feff diagrams are the best choice for 
such a comparison, because they rely on the properties which 
are most directly measured and are most easily interpreted. A 
good match to the properties of the dEB is found for an age 
of roughly 6Myr, a metal abundance ofZ^ 0.01, and a solar 
helium abundance. We have modelled the helium lines of both 
components of DW Car, using non-LTE calculations, to find a 
roughly solar or slightly subsolar helium abundance which is in 
acceptable agreement with that inferred from comparison with 
the theoretical models. The age is in good agreement with lit- 
erature values for Cr. 228, and the inferred metal and helium 
abundances are (to our knowledge) the first measurements for 
this cluster, illustrating the usefulness of detached eclipsing bi- 
naries for the investigation of stellar clusters. 
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8.1. DWCar, a Be star system 

The emission-line nature of the DW Car system is of interest 
both because of its implications for the properties of the system 
and in the study of classical Be stars. The Be nature of DW Car 
rests only on strong double-peaked emission in Ha and a possi- 
ble small emission in H/3: there is no noticable emission in the 
helium or metal lines of the system. The presence of double 
peaks indicates that the inclination of the circumbinary disc is 
high, which is consistent with the properties of the orbit of the 
eclipsing stars and the hypothesis that the orbit and disc should 
have a common orientation. 

The rotational velocities of the components of DW Car are 
towards the lower end of the distribution of Be star rotational 
velocities (e.g., Chauville et al. 2001 . Zorec, Fremat & Cidale 
12005 > . The significance of this is not straightforward to quan- 
tify, because we have not considered the effect of gravity dark- 
ening when measuring line widths (which tends cause the mea- 
sured rotational velocities to be underestimates), but is rein- 
forced because we have obtained an equatorial rotational ve- 
locity (which is the upper limit of the possible V sin / values). 

The Be phenomenon is thought to arise from the presence 
of some kind of pulsation in combination with a rotational ve- 
locity close to the critical limit for break-up of the star, causing 
the ejection of material to form a circumstellar disc (Porter & 
Rivinius'2003\ although magnetic fields may play a part. There 
are several effects which we would expect to see in DW Car. 
Firstly, about 90% of Be stars earUer than B5 exhibit small- 
scale variability (Porter & Rivinius"2003V There are no sugges- 
tions of this in our light curves of DW Car, but we would not 
have noticed anything with an amplitude significantly smaller 
than 5 mmag, which is the photometric precision of the v and b 
light curve data (Table |5J Paper 1). We would also expect that 
there would be some light contribution from the circumbinary 
material, but third light is negligible in the light curves. Both 
of these effects would have an insignificant effect on the light 
curve solution, so would not affect the reliability of the radius 
values. In the absence of intrinsic variability, the presence of a 
disc (which is confirmed by the Ha emission) must be caused 
by either dynamical effects due to the close-binary nature of 
DW Car or by capture of material from the surrounding envi- 
ronment. The unevolved status of DWCar is also consistent 
with the field population of high-mass Be stars (Zorec et al. 
120051 . 
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